Aharonov-Bohm-like effect for light propagating in nematics with 
disclinations 
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PACS 61 . 30 . Jf - Defects in liquid crystals 
PACS 42 . 70 . Df - Liquid crystals 

PACS 03 . 65 . Vf - Phases: geometric; dynamic or topological 

Abstract. - Using a geometric approach for the propagation of light in anisotropic media, we 
investigate what effect the director field of disclinations may have on the polarization state of 
light. Parallel transport around the defect, of the spinor describing the polarization, indicates the 
acquisition of a topological phase, in analogy with the Aharonov-Bohm effect. 
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Introduction. — In a recent article [1] two of us (CS 

and FM) studied the propagation of light in nematics 
with topological defects using a geometric approach due 
to Joets and Ribotta [2]. The model associates Fermat's 
principle of geometrical optics to geodesies in a hypothetic 
non-Euclidean space. This way, an effective geometry can 
be found which may be interpreted as "the cause" of the 
bending of the light paths, just as in general relativity. An 
important issue concerning this approach is the possibil- 
ity of using the mighty resources of differential geometry 
to calculate physical properties of such system. Indeed, in 
this Letter we show that spinors, which describe the polar- 
ization of a photon, when parallelly tranported around a 
loop encircling a disclination, acquire a topological phase 
, analogous to the Aharonov-Bohm effect [3] . 

EfTective geometry for disclinations in nematics. 

— We study wedge disclinations of generic strength k, 
whose director configurations, in the plane x — y, are given 

by [4] 

<p{e)=k0 + c, (1) 

where tp is the angle between the molecular axis and the 
a;-axis, 9 is the angular polar coordinate and c = f{0)- 
Selected director configurations can be seen on Figure 11.4 
of [4]. We assume the disclinations are straight and lie 
along the z-axis and the light rays propagate in the x-y 
plane so, effectively, we have a two-dimensional problem. 

Eq. (1) implies that f{2Tr) — (p{0) ~ 27rfc. Since the 
nematic molecules are symmetric under a rotation of n ra- 



dians around an axis perpendicular to the molecular axis, 
this result must be mr where n is a positive or negative 
integer. This means that k can only take integer or half- 
integer values. 

In ref. [1] we found that the effective geometry perceived 
by light traveling in the neighborhood of a disclination, in 
a nematic liquid crystal, is described by the line element 

+ {sin^ ^ + cos^ ^} p^d9^ 
- {2{a'^ - 1) sin cos ^} pdpdO, 



where ^ depends on the azimuthal angle 9 as 
^{9) = [k -l)9 + c 



(2) 



(3) 



and where a = ne/rio and rio and rie are the ordinary and 
extraordinary refractive indices, respectively. The radial 
coordinate r, used in ref. [1], has been rescaled to p — Uor. 

The above metric is associated to a curvature scalar 
given by 

2fc(fc- l)(l-a2) 



R = 



a^p^ 



■cos(20. 



(4) 



In fact, as we show below, this result is valid only for p 7^ 
since the idealized coreless defects will have a (5-function 
curvature singularity at p = 0. 
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Parallel transport of a vector in the effective ge- 
ometry. — Although our main interest in this work is 
the parallel transport of spinors, this section is important 
because it introduces the contribution of the singularity 

at the z-axis to the curvature, which justifies the main re- 
sult of this Letter. The parallel transport of a vector in 
a loop involving curved space is equivalent to a local ro- 
tation of the vector. The rotation angle contains precious 
information on the nature of the encircled curved space. 
We use this information to infer the contribution of the 
singularity at p = to the curvature. To do this we use 
the powerful and concise calculus of differential forms [5] . 

We introduce an appropriate dual 1-form basis 
(coframe), defined in terms of local zweibein fields by 
= e^^dx^, where 

= {cos^dp + sin ^ pdO) , (5a) 
= a {sin ^dp — cos ^pd6) . (5b) 

We introduce an affine spin connection 1-form and de- 
fine the torsion 2-form and the curvature 2-form, respec- 
tively, by 



where 



DO 

-n-fc — 



Tfe« A = de" + uj^A e^ (6a) 
Rbcde'^Ae''. (6b) 



A 
B 



a 



a 



-K sin^ cos^ , 



K a - 

- + K 

a a 



1 



■ cos ^. 



(11a) 
(lib) 



The spin connection Tg corresponds to orbits with the ra- 
dial coordinate p constant and Fp corresponds to curves 
with the 9 coordinate constant. The holonomy matrix 
U{j) associated with the parallel transport of vectors 
around a closed curve 7 is defined by 



U{'-f) = V exp 



dxT„ 



(12) 



where V is the ordering operator. The set of all holonomy 

matrices forms the holonomy group which is defined as 
the group of linear transformations of the tangent space 
TpM induced by parallel transport around loops based at 
a point p. This group possesses information regarding to 
the curvature of the effective geometry. 

We are interested in analyzing circular orbits. In this 
case the holonomy possess only the azimuthal contribution 
to the spin connection and is written as 



U{'ye) = "Pexp 



Tnd9 



(13) 



These are the Maurer-Cartan structure equations. Since 
we already know that the effective geometry is Rieman- 
nian [1] we make use of the torsion-free condition for 
this geometry, which implies that the first of the Maurer- 
Cartan equations becomes 



de" + uj'^Ae^ 



0. 



(7) 



From the 1-form basis we can determine the non-null con- 
nection 1-forms: 



sm t, cos ^ — 

a p 



■ cos^ C dO 



a 



(8) 



The connection 1-forms transform in the same way as 

the gauge potential of a non-Abelian gauge theory, which 
means that any two elements of the group do not com- 
mute. 

This result takes us to the following spin connection for 
the effective metric (2) 



and 



r. = 



B 
-B 



A 




(9) 



(10) 



which results in 



cos A — sin A 
sin A cos A 



(14) 



where 
A 



+ 



^ h QfTT^ k 

- 1 J sin [47r(A: 
a \ 



1) + 2cl - sin 2c 



4(fc - 1) 



k (15) 



The matrix U (7^) gives the effective rotation of the paral- 
lelly transported vector but it includes a 27r rotation due 
to the loop. Therefore x = 27r — A is the angle between 
the final and initial versions of the vector. 

Eq. (15) gets greatly simplified when the allowed values 
of k (integer or half-integer) are taken into account. For 
A; = 1 it is easily seen that 



X(fc = 1) = 27r(l - a) 
and for the remaining integer or half-integer k 



Xik 1) = 27r - TT 



1 

a + - 
a 



k. 



(16) 



(17) 



Notice that, when k = 1 and c = the metric (2) re- 
duces to that of a cosmological object, the cosmic string, 
as already reported in ref. [6] and, accordingly, the out- 
come of our calculations reduces to the know result for the 
cosmic string [7], where x = 27r(l — a). This, at first sight. 
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may seem incompatible with the null curvature obtained 
when one makes A; = 1 in equation (4). In fact, as it is 
well known in the cosmic string lore there is a (5-function 
curvature singularity at p = which is responsible for the 
change in the direction of the transported vector. As we 
show below this also happens for the generic k case. 

Before analysing the more general case we recall that 
for fc = 1 and c = the geometry is that of a cone [6]. 
Fig. 1 clearly shows the effect of the parallel transport of 
a vector alongside a circle in the Euclidean plane where a 
wedge of angle 27r(l — a) was removed in order to make 
the cone by identification of the lose edges. The resulting 
angle between the initial and final vectors is x = 27r(l — a), 
in agreement with eq. (16). 




Fig. 1: Parallel transport of a tangent vector around a circle 
concentric with the tip of a cone. The angle between the initial 
and final versions of the vector is x = 2tt{1 — a). 



Now, the above results have to be compatible with the 
Gauss-Bonnet theorem [8] which states that, for a regular 
Riemannian surface. 



KdA = 2tt 



Q 



dQ 



kg ds ^ 



(18) 



where K is the Gaussian curvature, dQ is the boundary 
of the singly-connected domain Q and kg is the geodesic 
curvature of the line dQ. Since our interest is circular 
orbits, we take Q to be a disc of radius p centered at 
p = Q and its circumference to be dQ. 

There is a very nice relationship between the Gauss- 
Bonnet theorem and the element of the holonomy group 
(14) giving the angle of local rotation of a vector after 
parallel transport in a loop [9] , 



X = 27r- 



dQ 



kgds = 



KdA, 



(19) 



which we explicitly verify below. 

For a curve in a surface with line element ds^ — Edv? 
2Fdudv + Gdv^ , the geodesic curvature is [8] 



kg = ^/EG - F2 [-Tj^ii^ + Tl^v^ - {2rj^ - Tl^)u'^i! 



(2rL 



T22)uv'^ + uv — vii] , 



(20) 



where Tjj are the Christoffel symbols of second kind and 
the derivatives are taken with respect to an affine param- 
eter along the curve, which we choose to be s itself. We 



choose u — 9 and v = p. The metric given by eq. (2) yields 

, we 

(21) 



%/ EG — = ap and, since dQ is a circle p — const, we 
have 



ds — y sin ^ + o? cos^ ^pdQ, 
from (2). This gives immediately 

Q = p"i(sin2 ^^o? cos^ ^Y^''^ 
The geodesic curvature reduces then to 

kg = ~apVl,B\ 
Now, using the metric (2) we obtain 



(22) 
(23) 



r fl 



Finally, 



P 



{a2(l_fc) + fc[l_ (i_a2)cos2^]2}. (24) 



kgds = 



It: 



(1- fc)a 



sin ^ -I- cos^ ^ 



with i{Q) = (fc 
to 



for k — 1 and 



fc 

-\sm'^^ + a'^cos'^n'>de, (25) 
a 



c. These integrals are easily solved 

(26) 



dQ 



kgds = 271 a 







1" 


/ kgds = TT 






'dQ 




a 



(27) 



for fc ^ 1 integer or half-integer. 

Since we are working in two-dimensional space, the 
curvature scalar is related to the Gaussian curvature by 
K — R/2. Using this, we substitute eq. (4) into the left 
hand side of eq. (18), with dA = apdpdO, which is easily 
obtained by taking the wedge product between and 
(eqs. (5a) and (5b)). For fc = 1 we have that K = R/2 = 
but X = 2T:a indicates that there should be a 5-function 
contribution to K. Therefore, in order that the Gauss- 
Bonnet theorem is satisfied. 



K = 



1 - a\ 5{p) 



a 



P 



(28) 



for fc = 1, in agreement with ref. [10]. 

Now, for fc ^ 1 the integral J KdA — again as it 
can be easily verified taking K = R/2 and using eq. (4). 
Again, it means that eq. (4) is incomplete since x = 27r — 
TT [a -I- i] fc. Therefore 



K 



1 


1 - 
















a 






a / 





fc(fc- l)(a2 - 1) 



2 9 



Sip) 



i{2[(fc- l)0 + c]} (29) 



for fc ^ 1 integer or half-integer. This means that, even 
though for p there is a p- and 6'-dependent curvature 
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associated to the disclination, the average of this curvature 
term in a disk centered in the origin is null. The singular- 
ity there is the sole responsible for the parallel transport 
effects. In real systems, the curvature (5-function singular- 
ity in the core of the defect will be smoothed out. But 
still, one can think of a localized flux of curvature ac- 
cross a very small area, the disclination core cross secion. 
So, even though the light ray does not crosses the core 
of the disclination, its spin angular momentum will "feel" 
the curvature flux, in a clear analogy with the Aharonov- 
Bohm effect. This justifies the next section on the parallel 
transport of a spinor. 

Parallel transport of a spinor in the effective ge- 
ometry. — An important question that emerges when 
we study the parallel transport of vectors is what happens 
with more complex fields when they also undergo a paral- 
lel transport. Vectors in a curved background have their 
orientation changed after a complete loop if there is a non- 
vanishing average curvature in the region surrounded by 
the loop. Spinors, as well, have their components changed 
if the same basis is used. This can be interpreted as a 
rotation in spinor space or the acquisition of a topological 
phase like Berry's [14]. 

Motivated by the spin angular momentum of the light, 
or polarization, we are interested in the study of the par- 
allel transport of a spinor in the effective geometry asso- 
ciated to the disclinations. The spinorial connection is 

rp(a:) = -^<7a7^ (30) 

where 7q are the flat-space Dirac matrices and lo"^ are the 
coefficients of the spin connection. The holonomy matrix 
is given, as in the parallel transport of vectors, by eq. (12) 
but in terms of the spinorial connection T^i{x) instead of 
the spin connection. 

For circular orbits we find the following expression for 
the spinorial connection in terms of the Dirac matrices 

^e{x) = -j^^le (717^ - 727^) 

= l^le^<^^ (31) 

where we have used the following representation for the 
7 matrices in terms of the Pauli matrices: 7^ = ia^ and 
7^ = —ia^. Again, the phase associated with the above 
spinorial connection is given by eq. (13). Integrating and 
expanding the exponential in eq. (13) we obtain 

C/(7e) = cos(x/2)/ + ^sin(x/2)fT^ (32) 

where / is the identity matrix, cr^ the Pauli matrix and x 
is given either by eq. (16) or eq. (17). Therefore, when a 
two-component spinor is parallelly transported in a closed 
orbit around the curvature flux given by eq. (28) or eq. 
(29) it undergoes a phase change given by 

l^') = C/(7.)|^)= ("'q^' J^/2)w^ (33) 



with phase angles corresponding, respectively, to eq.(16) 
and eq. (17). 

Discussion. — The phase (33) was obtained for a cir- 
cular path around the line defect. But this is a topo- 
logical effect and therefore the loop can be deformed 
without changing the result as far as the defect core is 
kept inside the loop. Moreover, since the line integral 
f — + Jb — Ia ^ Ia ^ following different paths (see 
Fig. 2), eq. (13) can be split into two parts, each inte- 
gral to be done on a separate path between points A and 
B. This way, the result (33) gives the polarization phase 



B 




A 



Fig. 2: Integration paths for eq. (13). 

difference between two light beams travelling in opposite 
sides of the defect. 

Considering a beam of linearly polarized light travelling 
past a disclination, if the input state is 

l^> = ^(l+> + |-», (34) 

the output state, obtained from eq. (33), is 

1^') = i= [exp(*x/2)|+) + exp(-zx/2)|-)] . (35) 

The internal product between these two states is 

m^')\^ = I [exp(jx/2) + exp(-»x/2)]' = cos2(x/2). 

(36) 

Then, according to Malus's law, the polarization plane 
of the light rotates by a relative angle x/2. That is, if 
the original beam is split in two, each of them travelling 
on either side of the defect, the relative rotation of the 
polarization plane is given by the angle x/2 when they are 
rejoined. Similar results apply to circularly and elliptically 
polarized beams, as can be easily verified. 

The rotation of the polarization in eq. (36) is analo- 
gous to what is observed when a linearly polarized laser 
beam propagates down a nonplanar optical fiber [15]. This 
effect has been interpreted as a manifestation of Berry's 
phase [19,20]. Like the original Aharonov-Bohm effect [3], 
where a flux of magnetic field provokes a phase shift in 
a split electron beam, the manifestation of Berry's phase 
can be associated to a flux of a gauge fleld in parameter 
space [21]. In our case, the topological defect is associated 
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to a flux of curvature in ordinary space. Since only the 
curvature associated to the defect core contributes, the 
phase is topological and should not depend on the local 
geometry. On the other hand, the defect is a particular 
configuration of nematic molecules. If one or a few of them 
have their orientation disturbed it is no longer the same 
defect. Then the phase associated to this deformed defect 
may not be the same as the one obtained here. Thermal 
agitation and molecular reorientation due to the electric 
field of a polarized laser beam are then effects that may 
affect experimental measurements of the predicted phase 
shift. To minimize thermal effects the sample should be as 
far below the isotropic to nematic transition temperature 
as possible. Since the reorientation angle is very small for 
laser beams with low power [11, 12] the minimum possi- 
ble power should be used. The closest the beam passes 
the core of the defect the better because the strong "an- 
choring" of the nematic molecules near the core will offer 
resistance to the torque provoked by the external field. 

In order to do an experiment to measure the effect pre- 
dicted here one needs disclinations in the nematic sample. 
These defects are spontaneously generated at the isotropic 
to nematic phase transition. In a recent article, Mukai and 
coworkers [13] produced disclinations in a lyotropic mix- 
ture (potassium laurate, decanol and water) in a sealed 
planar glass cell. Although their interest was in the for- 
mation statistics of the defects (which by the way repre- 
sents another analogy with cosmic strings) their procedure 
for generating and identifying the defects can be followed. 
The defects are observed in an optical microscope with 
the sample placed between crossed polarizers. If a binoc- 
ular microscope is used, one of the oculars can be used 
to focus the laser beam near a chosen defect and the sec- 
ond polarizer rotated in order to maximize the intensity 
at a photodiode conveniently placed after the second po- 
larizer. This will give the phase change in a single beam 
travelling past a disclination. One can try also an in- 
terference experiment with the aid of the converging lens 
behavior predicted for some of the disclinations [1] since 
they naturally force light rays following paths on opposing 
sides of the defect to meet. In this case the relative phase 
shift between the polarizations of the beams on different 
sides of the defect is measured. Although in the original 
Aharonov-Bohm effect there is no deflection of the beam, 
the analogy is in the fact that in both cases there is a 
confined field fiux causing a topological phase shift. As 
described above, the deflection is an advantage that can 
help in the detection of the effect. 

Diffraction of light by a disclination will certainly be af- 
fected by the effective geometry described here. Although 
the more interesting cases will be those of asymmetric de- 
fects (under present study) we can infer some properties of 
the k = I case (symmetric defect) from their cosmic ana- 
logues. In reference [6] it was shown that the k = 1 defect 
is analogous to a cosmic string in the sense that they share 
the same space geometry. Linct [16] and, more recently, 
Yamamoto and Tsunoda [17] investigated the propaga- 



tion of a plane wave past a cosmic string in the limit of 
geometrical optics. From their work we obtain that the 
cylindrical diffracted wave has an order of magnitude 
where A is the amplitude of the incident plane wave and 
a = For cosmic strings [18] a — 1 ~ iG^/c^ « 10~^, 
where G is the gravitational constant and /x is the mass 
density of the string, which for grand unification strings 
fi w lO^^g/cm. This indicates that, while the diffraction 
effect may be very difficult to observe for cosmic strings, 
it may be observable for disclinations in nematic liquid 
crystals which have a « 10*^. 

A recent experiment that is somehow related to the re- 
sults presented here was done by Marrucci and cowork- 
ers [22]. Injecting circularly polarized light in q plates, 
which are in fact cross sections of the defects discussed 
here, they observed the conversion of the spin angular 
momentum into orbital angular momentum. The inci- 
dent Gaussian beam is therefore converted to a helical 
mode preserving total angular momentum in this process. 
Their setup is exactly orthogonal to ours but nevertheless 
a change in the helicity is observed although not in the 
sense of our work. In our case there is only a rotation of 
the polarization without conversion of spin to angular mo- 
mentum therefore keeping a Gaussian beam in its original 
mode. 

A number of experiments have been carried out on the 
manifestation of Berry's phase for the photon polarization. 
Among them is the observation of the rotation of linearly 
polarized light propagating down a helically wound, single- 
mode optical fiber [15,20]. We have seen above that the 
polarization of light propagating in the neighborhood of 
disclinations is affected in a similar way by the director 
field of the defect. Using a geometric approach we showed 
that this polarization change is due to a flux of curva- 
ture associated to the core of the disclination. This is 
clearly a manifestation of Berry's phase and analogous to 
the Aharonov-Bohm effect in the sense that the spino- 
rial connection F^^ has the same role as the vector po- 
tential in provoking the phase change. Similar effects 
were predicted for disclinations in elastic solids [23] and 
for femionic quasiparticles in superfluids [24]. 
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